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ABSTRACT 

The total feedback energy injected into hot gas in galaxy clusters by central black holes can be 
estimated by comparing the potential energy of observed cluster gas profiles with the potential energy 
of non-radiating, feedback-free hot gas atmospheres resulting from gravitational collapse in clusters 
of the same total mass. Feedback energy from cluster-centered black holes expands the cluster gas, 
lowering the gas-to-dark matter mass ratio below the cosmic value. Feedback energy is unnecessarily 
delivered by radio-emitting jets to distant gas far beyond the cooling radius where the cooling time 
equals the cluster lifetime. For clusters of mass 4 — 11 x 10 14 M Q estimates of the total feedback 
energy, 1 — 3 x 10 63 ergs, far exceed feedback energies estimated from observations of X-ray cavities 
and shocks in the cluster gas, energies gained from supernovae, and energies lost from cluster gas by 
radiation. The time-averaged mean feedback luminosity is comparable to those of powerful quasars, 
implying that some significant fraction of this energy may arise from the spin of the black hole. The 
universal entropy profile in feedback-free gaseous atmospheres in NFW cluster halos can be recovered 
by multiplying the observed gas entropy profile of any relaxed cluster by a factor involving the gas 
fraction profile. While the feedback energy and associated mass outflow in the clusters we consider 
far exceed that necessary to stop cooling inflow, the time-averaged mass outflow at the cooling radius 
almost exactly balances the mass that cools within this radius, an essential condition to shut down 
cluster cooling flows. 

Subject headings: X-rays: galaxies: clusters; black hole physics 



1. INTRODUCTION 

Massive black holes in the cores of group and cluster- 
centered galaxies become active when a tiny mass of clus- 
ter gas is accreted, causing large amounts of "feedback" 
energy to be deposited in the surrounding cluster gas. 
The energy provided by occasional feedback events, in 
the form of jets containing supra-thermal and relativis- 
tic particles, increases the entropy of both nearby and 
very distant cluster gas. Feedback jets create expanding, 
shock-driving cavities in the hot cluster gas, increasing 
its entropy and delaying its flow toward the black hole by 
radiation losses. In the absence of feedback, the loss of 
entropy by thermal X-radiation is accompanied by a slow 
cooling inflow toward the black hole which, if unabated 
over time, would concentrate a mass of centrally cooled 
gas in or near the central black hole far exceeding limits 
set by stellar velocities and other observations. For ex- 
ample, intermittent cavities formed from feedback events 
of energy 10 59 ergs every 2 x 10 8 years at 10 kpc can arrest 
the currently observed gas density and temperature pro- 
files (and low central cooling) in the Virgo cluster for ~ 3 
Gyrs (Mathews 2009). Including more distant feedback 
events at 50 kpc can maintain the observed gas profiles 
and low cooling rate for several additional Gyrs. Entropy 
and cosmic rays delivered to the cluster gas by feedback 
episodes cause gas to flow out in the cluster potential, 
offsetting the cooling inflow. 

In general, however, central black holes are unable to 

University of California Observatories/Lick Observatory, De- 
partment of Astronomy and Astrophysics, University of Califor- 
nia, Santa Cruz, CA 95064 mathews@ucolick.org 



provide the optimal or minimum feedback energy that 
must be deposited at every radius in the cluster gas just 
sufficient to shut down black hole accretion of locally 
cooling gas. Instead, black holes typically over-react in 
a clumsy fashion, depositing much more energy than the 
minimum required, much of it in distant regions of the 
cluster where the radiative cooling time exceeds the age 
of the cluster. After a few 10 8 yrs following a feedback 
event most of the feedback energy converts to potential 
energy, as the entire gaseous cluster atmosphere adia- 
batically expands outward (Mathews & Brighenti 2008). 
Consequently, the increased cluster gas entropy is nec- 
essarily related to a reduction of the cluster gas density 
relative to the local dark matter. 

We discuss here an approximate estimation of the total 
feedback energy received by cluster gas during the clus- 
ter lifetime by comparing gas potential energy profiles in 
observed clusters with that of idealized gas density distri- 
butions resulting from "adiabatic" gravitational collapse 
into the cluster halo in the absence of radiative cool- 
ing and associated feedback. We show that this energy, 
~ 10 63 ergs, far exceeds the energy lost by radiation dur- 
ing the cluster lifetime and consequently the minimum 
energy required merely to stop the cooling inflow. The 
collective energy from all supernovae also provides a neg- 
ligible fraction of the total feedback energy. Since only 
10 percent of the hot baryonic gas in massive clusters 
cools to form stars, star formation can also be ignored 
in our estimate of the global energetics of the cluster gas 
where we seek an overall accuracy of ~ 25%. 

Estimates of the total feedback energy from cluster- 
centered black holes are possible for two reasons. First, 
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the potential energy of cluster gas in hydrostatic equi- 
librium in a fixed dark halo potential can be found by 
integrating outward from the cluster center. Second, to 
a good approximation, the formation of dark halos and 
their gravitational potential proceeds from the inside out- 
ward as the size of the virialized region in the dark halo 
increases with time. In view of this latter point, it is 
possible to estimate the increase in potential energy of 
the cluster gas due to feedback without knowing when 
the feedback occurred. 

We begin by estimating the current gas density, tem- 
perature, mass, entropy and potential energy profiles ex- 
pected at zero redshift in the absence of radiative losses, 
star formation, and black hole feedback, referred to as 
the "adiabatic" cluster atmosphere. Then, assuming hy- 
drostatic equilibrium, potential energy profiles are esti- 
mated from gas distributions in observed clusters having 
the same total mass. When comparing the difference be- 
tween the potential energy evaluated at the same mass of 
cluster gas with and without feedback, we find that hot 
gas in the idealized adiabatic atmosphere must expand 
significantly to resemble cluster gas profiles currently ob- 
served. The total feedback energy associated with this 
expansion, ~ 1 — 3 x 10 63 ergs, comfortably exceeds en- 
ergies of the most powerful known individual feedback 
events. Central to our feedback energy estimate is the in- 
verse relationship between increasing cluster gas entropy 
and decreasing gas fraction, the ratio of gas to total clus- 
ter densities, which is lowered by a global expansion of 
the cluster gas. 

If this feedback energy is released during periods of 
radiatively efficient central accretion with ~ 0.1 of the 
accreted mass returned as feedback energy to the clus- 
ter gas, we find that the final black hole masses in large 
clusters would exceed those observed. Alternatively, the 
feedback energy may be provided from the rotational en- 
ergy of rapidly spinning central black holes. 

2. ESTIMATE OF TOTAL FEEDBACK ENERGY 

A key element in our estimate is an assertion that adi- 
abatic (non-radiating, non-feedback) gaseous cluster at- 
mospheres can be fit with the same properly normalized 
NFW profile as the dark matter. Among the many com- 
putations of galaxy cluster formation that include both 
baryons and dark matter, surprisingly few have analyzed 
in detail the deep similarity between the final adiabatic 
gas structure and the radial NFW distribution of dark 
matter. Notable exceptions to this are the density and 
entropy profiles in dark halos and adiabatic gas described 
in the cosmological cluster formation calculations of Eke, 
Navarro & Frenk (1998) and in particular Faltenbacher 
et al. (2007). The cluster gas entropy can be character- 
ized with S g = cr 2 //? 2 / 3 where p is the gas density and 
a = [3fcT/(/zm p )] 1/ ' 2 is the thermal velocity dispersion. 
By analogy, Faltenbacher et al. define a corresponding 
dark matter entropy Sdm = cr 2 /p 2 ^ 3 for which p is the 
dark matter density and a is the 3D velocity dispersion of 
collisionless dark matter particles. The detailed cosmo- 
logical cluster calculations of Faltenbacher et al. (2007) 
using GADGET2 reveal that, apart from a small central 
gaseous core, the gas density and entropy profiles in the 
adiabatic case share identical, appropriately scaled NFW 
profiles regardless of cluster mass. Beyond the central 



core where S g > Sdm, both S g and Sdm vary as power- 
laws, cx r , but with slightly different normalizations, 
Sg/Sdm = 0.71 ± 0.18. However, as discussed by Fal- 
tenbacher et al., the final zero redshift gas configuration 
in the gas contains residual subsonic macroscopic veloci- 
ties which, if damped, would bring the gas and dark mat- 
ter entropies into near perfect agreement, S g /Sd m ~ 1- 
We adopt the assumption, implicit in the discussion of 
Faltenbacher et al., that the two entropies are in fact 
equal but the gas remains undamped because the ap- 
propriate physical damping mechanisms are absent from 
GADGET2. The artificial viscosity in this code that 
damps the accretion shock may require an unphysically 
long time to damp the residual subsonic velocity field. 
Assuming S g ~ Sdm is also more consistent with the 
smaller ~ 5% component of turbulent energy found for 
relaxed clusters in recent cosmological simulations using 
the ENZO code (Vazza et al. 2011). Finally, cluster virial 
masses are determined from observed gas pressure pro- 
files without allowing for undamped kinetic energy and 
its associated pressure. This assumption (which underes- 
timates the virial mass) is also consistent with S g « Sdm- 
While it is comforting that the radial profiles of adia- 
batic gas and dark matter are nearly identical, this may 
nonetheless be surprising due to the very different na- 
ture of their dissipative mechanising. Of more relevance 
to our discussion is the similar NFW shape of adiabatic 
gas and dark matter density profiles. The assumption 
S g ~ Sdm ensures that the adiabatic gas has experienced 
the same dissipative history as the dark matter. 

The radial dark matter NFW distribution is shaped 
by all entropy-producing dissipations that occurred dur- 
ing both smooth accretion of diffuse dark matter as well 
as inhomogeneous accretion of smaller groups, clusters 
and filaments that merged into the final cluster poten- 
tial. The entropy increases that accompanied the forma- 
tion of these smaller structures are also embedded in the 
final NFW structure. Furthermore, the results of Fal- 
tenbacher et al. (2007) indicate that all dissipative infor- 
mation about the cluster merger history is also encoded 
in the NFW gas density profiles in idealized adiabatic 
baryon atmospheres formed by gravity alone. 

When computing purely gravitational (adiabatic) col- 
lapse in cosmological calculations without radiation or 
feedback, a controversy has arisen in recent years con- 
cerning differences in the dissipation between dark mat- 
ter and baryons in cluster cores where S g > Sdm (recently 
reviewed by Borgani & Kravtsov 2009; Springel 2010; 
Vazza 2011). In mesh-based calculations the baryons are 
found to have large, well-resolved central density cores 
which are larger than those computed with smoothed 
particle hydrodynamics. Possible origins for this discrep- 
ancy have been proposed and discussed by Mitchell et al. 

2 In adiabatic cluster simulations gases of different entropies ev- 
idently mix in the cluster core, raising the total gas entropy, par- 
ticularly in grid-based computations. However, beyond the core 
we expect the density of gas and dark matter to share the same 
appropriately normalized NFW profile. This profile similarity has 
been verified in many calculations including the Santa Barbara 
cluster, an average of 12 different structure formation codes using 
identical cosmologies and initial conditions (Frenk et al. 1999). In 
more recent adiabatic simulations the effective dark matter temper- 
ature profile Td m { T ) = (^'m p /3k)a 2 has been shown to be identical 
(within 10%) to the gas temperature profile T(r) (e.g. Host et al. 
2009; ZuHone 2011). 
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(2009) and Vazza (2011). For our purposes here we as- 
sume two limiting adiabatic gas density profiles following 
pure gravitational baryonic collapse into dark halos: core 
and no core. 

2.1. Adiabatic cluster gas atmospheres without cores 

Consider first the "no core" case in which baryons 
and dark matter suffer the same dissipation so the post- 
collapse gas density profile is identical to that of the total 
density but scaled down by the universal baryon fraction 
fb = 6.17, i.e. p = f b p t (e.g. Faltenbacher et al. 2007) 
where p t is the total density dominated by dark mat- 
ter. We seek the idealized radial structure of cluster gas 
formed without feedback or radiative losses and which 
has evolved to the current time. The temperature profile 
can be found by integrating the equation of hydrostatic 
equilibrium 

d0 J ldp\ 

9 (1) 



dr 



p dr 



where 6 = kT/pm p and g = GM{r)/r 2 . The self- 
gravitation of the gas is implicitly included in the total 
mass M(r) that includes both dark matter and gas. The 
total matter density and mass profile in a cluster with 
virial mass M v are given by the usual NFW relations, 
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and M(r) = M v 
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y(i + y) 2 w /(c) 

where c = y/(r/r v ) is the concentration and 

f(y) = \n(l + y)-y/(l + y) S c = A c c 3 /3/(c). (3) 

The radius where the local cluster density is A times 
larger than the critical density p c = 3H 2 /8nG — 9.24 x 

10 -30 g cm -3 ; s 



r A = {3M v /4nAp c ) 



1/3 
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and the virial radius is r v = r\ c where A c = 
178(fi m ) 45 = 103 (Eke, Navarro, Frenk 1998) when 
Ojy = 0.3. For nearby clusters we assume H = 70 km s _1 
Mpc -1 . Finally, we adopt a convenient concentration- 
M v relation c = 9(M„/lO 14 M )-°- 172 for low redshift 
clusters (Buote et al. 2007). 

In our integrations of equation (1) using the observed 
gas density we consider two clusters each of which are 
composites of two very similar clusters selected from 
the sample of nearby relaxed clusters observed with the 
Chandra telescope by Vikhlinin et al. (2006). For our 
purposes we choose massive clusters in which the full 
impact of feedback energy can be absorbed and which 
do not presently have huge cavities that dominate and 
disrupt the cluster gas profiles observed. Of particular 
interest are the NFW parameters determined for these 
clusters and values of the observed cluster gas fraction 
f g = pj pt plotted in Figures 3-14 of Vikhlinin et al. To 
avoid being distracted by spurious observational errors 
and to improve the quality of our estimates, we com- 
bine the mean properties of two pairs of clusters that 
share nearly the same virial mass: (A133 + A383) and 
(A478 + A1413). The two averaged clusters are subse- 
quently referred to as composite clusters 1 and 2 respec- 
tively. Relevant observed NFW properties of all four 
clusters and the the two composite clusters are listed in 



Table 1. Each composite cluster is chosen to have values 
of r 50 o and M 50 o that are averages of the two observed 
clusters. Composite clusters 1 and 2 have virial masses 
M v = 4.34 x 10 14 and M v = 11.3 x 10 14 M Q respectively. 

Our first integrations of equation (1) are for adiabatic 
"no core" versions of the two composite clusters using 
p( r ) = fbpt(r) for the gas density. We choose an initial 
temperature at some very small cluster radius and seek 
solutions for which the gas entropy S = 9/p 2 / 3 varies like 
a self-similar power law at large radius S oc r , simi- 
lar to entropy profiles in detailed adiabatic gravitational 
collapse computations (e.g. Tozzi & Norman 2001; Voit 
et al. 2005). This is a well-defined initial value problem 
with unique solutions. By varying the initial tempera- 
ture, an exact value and solution can be found for which 
S oc r q where q « 1.2 remains constant over a large 
range of cluster radius near the virial radius and beyond. 
The value of q emerges naturally from the solution as an 
eigenvalue that is not imposed in advance. 

Solid lines in Figures 1 and 2 show adiabatic post- 
collapse profiles of gas temperature, entropy and den- 
sity respectively for composite clusters 1 and 2 found by 
integrating equation (1). The hydrostatic post-collapse 
cluster gas atmospheres in Figures 1 and 2 have a broad 
temperature maximum near r/r v =0.1 and distant en- 
tropy profiles S oc r 1,3 that are very similar to the Tozzi- 
Norman computation. Our simple approximation closely 
resembles sophisticated cluster gas profiles computed in 
detailed cosmological simulations that include mergers 
of smaller bound systems and gas inflowing in filaments 
(e.g. Borgani & Kravtsov 2009). Gas and dark matter 
in these merging sub-systems and filaments have already 
experienced some dissipation when they enter the cluster 
virial radius. Apart from this we consider no other type 
of hypothetical ad hoc pre-heating. 

2.2. Observed cluster gas atmospheres without cores 

Next we perform similar integrations using the ob- 
served cluster gas density profiles p a b s (r) , approximately 
including the gas self-gravity in the total NFW mass dis- 
tribution. The data points shown with open circles in 
Figures 1 and 2 show observed gas densities found from 
Pobs{r) = fg(r)pt(r) for both clusters in each composite 
cluster at radii for which f g is determined by Vikhlinin 
et al. In addition we show (with filled circles) several ad- 
ditional approximate values measured directly from the 
Pobs(r) profiles plotted by Vikhlinin et al. (Closed and 
open circles have sizes roughly comparable to observa- 
tional errors.) 

The mean observed gas density profiles of our two com- 
posite clusters are fit with an analytic curve 



Pobs(r) 



fbScPc 



y 1 - a (yo + y) 2+a 



(5) 



where y = c(r/r v ). By design, p bs{r) asymptotically ap- 
proaches fbPt(r) = 0.17p t (r) as the radius continues be- 
yond the observed region, i.e. r > 0.5r v . As the cluster 
gas conserves baryons during feedback expansion, regions 
OI fg < fb in the cluster gas observed within ~ 0.5r„ 
must be compensated by regions of f g > fb in more dis- 
tant cluster gas. However, due to the increased volume 
available in the outer regions of the clusters, we expect 
the excess f g — fb to be much less than and difficult 
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to observe. Profiles for p bs(r) are shown with dotted 
lines in the bottom panels of Figures 1-4 where we take 
2/0 = 1-8 and 1.5 for cluster 1 and 2 respectively and 
a = 0.05 for both composite clusters. 

For the observed hot gas atmospheres in clusters 1 and 
2 equation (1) is solved in the same manner as before 
but with d\np bs/dr on the right hand side. The cor- 
responding hydrostatic temperature T(r) and entropy 
S(r) profiles are shown with dotted lines in the upper 
and central panels of Figures 1 and 2 for each compos- 
ite cluster. Both the temperature maximum and the 
somewhat flatter entropy profile that asymptotically ap- 
proaches the adiabatic profile resemble typical cluster 
observations (e.g. Vikhlinin et al. 2006; Pratt et al. 
2010). The total bolometric X-ray luminosity within the 
virial radius L x (r v ) = (3.83 x 10 44 , 1.58 x 10 45 ) erg s" 1 
respectively for composite clusters (1,2) are consistent 
with the L x — M scaling relation for observed clusters 
(e.g. Maughan 2007). 

While the two cluster atmospheres plotted in Figures 
1 and 2 - adiabatic and observed - are final zero-redshift 
gas distributions, we can also regard them as initial and 
final configurations of the cluster gas before and after 
feedback energy is deposited. Of particular interest is 
the potential energy of the gas at various radii in the 
two atmospheres that enclose the same integrated gas 
mass M g (r) . The potential energy of gas in the cluster 
is r 

PE(r) = / (f>(r) P 4irr 2 dr (6) 
Jo 

where the NFW potential 

GMy ln(l + y) 

^ rWT~ (7) 

also satisfies g = —d<j>/dr. In evaluating the potential 
<p(r) we again approximate the self-gravity of the gas by 
using the total cluster mass, and we ignore the small 
adiabatic expansion of the dark matter as gas moves out 
in the cluster potential. The NFW potential is expected 
to be valid within about twice the virial radius (Cuesta 
et al. 2008; Tavio et al. 2008). 

Note that the potential energy profile is determined by 
an integration out from the cluster center. Furthermore, 
to a good approximation, the cluster dark matter density 
and potential profiles remain fixed with time as the the 
size and mass of the virialized region increases with time, 
i.e. the mass distribution in virialized dark halos forms 
from the inside out. The inside-out character of halo 
formation is apparent for example from Figure 1 of Die- 
mand, Kuhlen, Madau (2007) who plot the accumulated 
density within many Lagrangian mass zones with time as 
a dark halo grows in size and mass. In view of this, our 
estimate of the change in PE due to non-gravitational 
feedback, assuming an unchanging total NFW mass dis- 
tribution, is independent of the time when the feedback 
occurred. As the dark halo grows, the global outflow of 
cluster gas driven by feedback energy can relocate the 
accretion shock beyond the instantaneous virial radius. 
But we do not consider extremely distant feedback events 
(r > 1.5 — 2r v ) that energize gas before it reaches the 
accretion shock and where the NFW potential may no 
longer apply. This restriction on feedback seems rea- 
sonable, particularly since the virial radius continuously 



increases with time and it is consistent with the gradual 
radial increase in f g toward fb observed in cluster gas. 

Radii and potential energy PE values at three loca- 
tions in the "no core" adiabatic and observed solutions 
having the same M g (r) are listed at the top of Table 2 for 
clusters 1 and 2. It is seen that the radius that encloses 
all the gas within the virial radius in the pre-feedback 
adiabatic cluster increases by about a factor 1.9 (clus- 
ter 1) or 1.5 (cluster 2) as the cluster gas expands due 
to feedback energy. This expansion has been verified by 
George et al. (2009) who observe a galaxy cluster with 
virialized, X-ray emitting gas well beyond the virial ra- 
dius. The increase in potential energy due to feedback 
when integrated to the same gas mass, shown as APE 
in Table 2, is huge, 1 — 3 x 10 63 ergs, greatly exceed- 
ing the most energetic known individual bipolar feed- 
back events < 10 62 ergs (McNamara et al. 2005; Guo 
& Mathews 2010a), which are sufficient to convert cool- 
core to non-cool-core clusters (Guo & Mathews 2010b). 
The mean luminosity of this enormous feedback energy, 
if spread over a typical cluster lifetime t c i ~ 7 Gyrs, is 
Lfb « 4 x 10 45 (|RE|/10 63 ergs) erg s _1 , comparable to 
continuously active quasars. 

2.3. Adiabatic and observed cluster gas atmospheres 
with cores 

In a second series of similar calculations, we consider 
adiabatic collapse atmospheres in which the baryons ex- 
perience an additional central dissipation that produces a 
density core. For this purpose, we adopt the recent mesh- 
based purely gravitational cluster collapse computations 
of Vazza (2011). We create a core from the no core den- 
sity profile by flattening the central gas density so that 
d\ogp/d\ogr « —2 at r/r v « 1.4 as in Figure 3 of Vazza 
(2011). However, to satisfy the stability requirement of 
radially increasing entropy, the gas density in the core 
cannot be perfectly flat. We find positive dS/dr when 
the core gas density gradually decreases as p oc r~ 0,3 and 
this solution is shown as a solid line in the bottom panels 
of Figures 3 and 4. In our approximation baryons in the 
core region are simply removed without a corresponding 
increase in f g beyond the core region. Corrections of 
this sort that must occur for baryon conservation, and 
which arc not apparent in the mesh-based calculations, 
are not important for our estimate here since the baryon 
fraction in the cored adiabatic atmosphere is essentially 
fb = 0.17 at large cluster radii. Finally, we assume that 
adiabatic cored density profiles are a universal function, 
scaling with r/r v for clusters of different virial masses. 
Performing similar integrations of the hydrostatic equa- 
tion (1), temperature and entropy profiles corresponding 
to the cored density profile are shown as solid lines for 
composite clusters 1 and 2 in the upper two panels of 
Figures 3 and 4 where it is seen that the entropy has a 
slowly sloping entropy "floor" , as expected. 

Repeating the same procedure as before, three pairs 
of equal gas mass M g {r) locations in the two cluster gas 
atmospheres - adiabatic "core" collapse and observed - 
are listed at the bottom of Table 2. When baryonic cores 
are present the potential energy increases APE are only 
slightly less than if the core is completely ignored. 

2.4. Further observational implications 
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Both of the starless adiabatic cluster atmospheres we 
consider - with and without cores - have local gas frac- 
tions that agree with the cosmic baryon value fb — 
0.17 near r v . However, the large expansion of cluster 
gas caused by feedback reduces the local gas fraction 
throughout the observed post-feedback cluster. The ob- 
served gas fraction at the most distant density obser- 
vations near r « 0Ar v is f g — 0.082, only about half 
the cosmic value. This shortfall is typical of all mas- 
sive clusters. The average value of the gas fraction 
at rsoo ~ 0.5r„ for clusters having mean temperatures 
greater than 4 keV is / s (r 500 ) ~ 0.12 ± 0.02 (McCarthy, 
Bower & Balogh 2007). When feedback-energized clus- 
ter gas expands out beyond r 5no , baryon conservation 
requires that f g > fb somewhere beyond r 500 , although 
baryon excesses have not yet been observed. Neverthe- 
less, it is possible that nearly baryonically closed systems 
in which APE << \PE(r v ) \ and f g (r v ) « 0.17 may exist 
in some fossil groups where the total mass of gas cooling 
near the central black hole is much less than in the more 
massive clusters considered here (Mathews et al. 2005). 

The final pair of entries for each of the four clus- 
ter atmospheres in Table 2 are designed for the major- 
ity of current cluster observations that extend only to 
~ ^500 ~ 0.57\,. Values of APE(r ac i) with r a d ~ 0.3r„ in 
Table 2 represent the minimum feedback energy consis- 
tent with clusters observed to ~ 0.5?v This minimum 
feedback energy - about 3 — 13 x 10 62 ergs from Table 2 - 
would be correct if f g increases abruptly to fb = 0.17 just 
beyond the outermost gas density observed. While a few 
recent cluster X-ray observations with the Suzaku tele- 
scope extend to the virial radius and beyond (e.g. George 
et al. 2009), these data may be more difficult to inter- 
pret. Detailed calculations of the cosmological evolution 
of massive clusters often indicate significant undamped 
subsonic gas velocities beyond about rsoo ~ 0.5r„, which 
may degrade the assumption of hydrostatic equilibrium 
(Evrard et al. 1996). Baryonic clumps or compressions 
occur in this same outer region, causing the observed 
gas density to appear too high and the entropy too low 
(Nagai et al. 2007a,b). 

2.5. Influence of stellar baryons 

In our estimates of the feedback energy we have ig- 
nored the small fraction of baryons in our clusters that 
convert to stars. Nevertheless, it is interesting to com- 
pare the mass of stellar baryons to the total baryonic gas 
mass that flows out at large cluster radii due to feed- 
back energy. The mean cluster stellar mass fraction /* 
decreases with cluster mass as /» oc M -0 - 55 (Andreon 
2010) but this variation has considerable cosmic scatter. 
For our clusters 1 and 2 with mass logM2oo = 14.56 
and 14,96 in M , the stellar fraction, /* w 0.008 - 0.015 
(from Figure 7 of Andreon 2010), is dominated by cosmic 
scatter. Adopting a mean value for this range of cluster 
mass, /» w 0.010, the total stellar mass within 0.5r„ in 
clusters (1,2) is (0.26,0.65) x 10 13 M©. By comparison, 
from Table 2 the total mass flowing out beyond 0.5r„ for 
(core, no core) versions of cluster 1 is (1.823, 1.983) x 10 13 
M Q and (3.69, 3.91) x 10 13 M Q for cluster 2. The frac- 
tion of outflowing gas mass consumed by star formation, 
^0.13 — 0.17, is another small correction that we neglect 
here. 



2.6. Specific feedback energy 

Our values of APE(r v ) = (8.34,29.1) x 10 62 ergs for 
cored clusters (1,2) translate into mean specific feedback 
energies APE/(0A7M V ) = (5.65,7.57) x 10 15 erg gm -1 
when applied to the total gas mass within r v in the orig- 
inal adiabatic clusters that extends out to (1.5, 1.8)r v in 
the post-feedback atmosphere. The specific feedback en- 
ergy can also be expressed as (3.6,4.8) keV particle -1 . 
These values appear comparable to those estimated from 
the Lx — T scaling relation for clusters by Wu, Fabian 
& Nulscn (2000), - 1 - 3 keV particle -1 . However, Wu 
et al. only consider particles in the mass of observed 
gas within r2oo! our § as masses are larger by (2.0,1.8) so 
our estimated feedback energies are more than twice as 
large as those of Wu et al. By comparison, our specific 
feedback energies are considerably less than those of Mc- 
Carthy, Bower, & Balogh (2007), > 10 keV particle -1 , 
probably because they use an extrapolation procedure 
for the observed gas density p bs(r) quite different from 
our equation (5). 

2.7. Energetics of radiation loss and supernovae 

We assume that the energy and entropy lost by radi- 
ation can be ignored. To verify this, consider the total 
energy radiated by the two composite clusters within the 
cooling radius E ra d ~ L x (r coo i)t c i during the cluster life- 
time t c i — 7 Gyrs. For composite clusters (1,2) we find 
rcooi = (97.7,119.8) kpc, L x {r cool ) = (1.36,4.96) x 10 44 
erg s -1 , and E rad = (0.30, 1.1) x 10 62 ergs. The total ra- 
diated energy for both clusters is a small fraction of the 
change in potential energy APE listed in Table 2, partic- 
ularly those values of the total APE(r v ) evaluated at the 
virial radius in the adiabatic solutions. Radiation losses, 
while not a major factor in the overall energy budget for 
our clusters, are nevertheless essential in driving the cen- 
tral cooling accretion that creates the feedback energy. 
The cooling radius (\og(r coo i/r v ) w —1.3) is small com- 
pared to the more extended regions in the bottom panels 
of Figures 1-4 where the cluster density is observed to 
be significantly below the maximum cosmic baryon den- 
sity, p b s < fbPt- Feedback energy from cluster-centered 
black holes has obviously caused huge outflows, remov- 
ing cluster gas from regions far beyond the cluster cooling 
radius. 

Feedback from Type II and la supernovae can also be 
ignored. For a single-population Salpcter IMF, about 
0.007 Type II supernovae events occur per solar mass of 
gas formed into stars. Consider the total mass of gas ob- 
served in the composite observed clusters (1,2), extrap- 
olated to the virial radius, about M g (r v ) = (4.7, 13.9) x 
10 13 M Q . If 10% of the baryons in this mass forms into 
stars, the total feedback energy from SnII, each of en- 
ergy 10 51 ergs, is of order Es n ii ~ (0.3,0.9) x 10 62 ergs 
which is substantially less than the APE(r v ) in Table 
2 evaluated at radius r v in the pre-feedback gas. If the 
average iron abundance in cluster gas is about 0.3 solar, 
the total mass of iron within r v in the clusters is roughly 
0.3 x 0.0017 xM g (r v ) = (2.4, 7.1) x 10 10 M where 0.0017 
is the fraction of iron by mass in the solar photosphere. 
If all the iron is created in Type la supernovae, each hav- 
ing energy 10 51 ergs and providing 0.7 M in iron, the 
total feedback energy from Type la supernovae cannot 
exceed about Es n ia ~ (0.3,1.0) x 10 62 ergs, which are 
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also less than APE(r v ) in Table 2. 

While AGN feedback dominates cluster energetics, su- 
pernovae may nevertheless contribute 10-20% of the to- 
tal feedback energy. Nagai et al. (2007a, b) computed a 
variety of gaseous atmospheres in clusters including su- 
pernovae of all types but without AGN feedback, as com- 
monly assumed. In their clusters the baryon fraction and 
entropy profiles are in reasonably good agreement with 
those observed by Vikhlinin et al. (2006) outside the 
central region, r > 0.2r/r v . But central overcooling in 
r ^ 0.2r/r v is a serious problem. At zero redshift in the 
models of Nagai et al. about 40% of the baryons within 
r ~ 0.5r„ are in the form of a centrally concentrated mass 
of stars and cold gas. The mass of this central concentra- 
tion of radiatively cooled baryons causes the gas density 
and temperature to peak up near the cluster centers un- 
like the observations. A large amount of AGN feedback 
energy, similar to that estimated here, is essential to re- 
move this overcooling gas before it forms into stars and 
relocate it to distant regions of the cluster. The outward 
flow of cluster gas that results from the creation of X-ray 
cavities is described by Mathews & Brighenti (2008) and 
Mathews (2009). AGN feedback has just begun to be in- 
cluded in recent cosmological cluster calculations where 
the overcooling problem is greatly alleviated (Teyssier et 
al. 2010; Puchwein et al. 2010; McCarthy et al. 2010). 

2.8. Feedback stops cooling flows within the cooling 
radius 

The approximate rate that mass cools in the two com- 
posite clusters in the absence of feedback can be esti- 
mated from the bolometric X-ray luminosity at the cool- 
ing radius, 

L x {r coo i) = M cf - = 1.50 xlO 41 — ^ {- T keV 

2pm p y M Q yr- 1 J 

(8) 

(in erg s ) where the cooling radius r coo ; is defined as 
that radius at which the local gas cooling time is equal 
to a typical cluster age t c i ~ 7 Gyrs. For cluster (1,2) 
with r cool = (98,120) kpc and L x {r cool ) = (1.36,4.96) x 

10 44 erg s _1 (determined with p bs(r)) we find M c f « 
(225,465) Mq yr -1 . This is the approximate rate that 
the cluster gas would cool in the absence of feedback 
energy. The ratio for the two clusters M c f t \/M c f t 2 = 
0.48 is comparable to but larger than the ratio of total 
feedback energies, APE\/ APE 2 ~ 0.29 for gas initially 
at r v in the "no core" solutions. 

As shown in Mathews (2009), feedback energy can 
create cluster outflows that balance and arrest cooling 
inflows caused by radiation losses. For this to hap- 
pen within a strongly radiating cluster core, the time- 
averaged rate that cluster gas flows out due to feedback 
expansion must be equal to the average mass inflow rate 
M c j due to radiative losses. First imagine the expanding 
flow of cluster gas due to feedback, ignoring energy losses 
due to radiation. In this limiting case the mass of gas 
within any radius r decreases because of feedback expan- 
sion from M a d{r) to M b(r) as the initial adiabatic at- 
mosphere evolves to the observed one during the cluster 
lifetime t c \. The mean positive feedback mass flow past 
radius r, AM(r)/t c i = [M ad {r) - M ohs {r)]/t cU is plotted 



in Figure 5 for both composite clusters. At the other 
limit when radiative losses are present but feedback is 
absent, the rate that gas cools within radius r, M c f(r) is 
related by equation (8) to the total rate that energy is ra- 
diated within this radius, L x (r). In particular, the mass 
of gas within the cooling radius in the observed cluster 
atmosphere, by definition of r coo i, is expected to flow in- 
ward at an average rate M c f(r coo i). Consequently, in the 
presence of both radiative losses and feedback gains, feed- 
back from the central black hole is expected to increase 
until the time-averaged rate that gas mass is driven out- 
ward past r coo i balances the average rate that mass would 
cool and flow inward past r coo i over the cluster lifetime. 

This balance is shown in Figure 5 where the total time- 
averaged rate that mass flows out past the cooling ra- 
dius, AM(r coo i)/t c i very nearly matches the approximate 
mean inflow rate M c f(r coo i) expected (from eqn. 8) at 
radius r coo i during time t c i if radiative cooling within this 
radius is uninhibited by feedback. Since inflow M c f(r) 
(estimated from L x (r) without feedback) and outflow 
AM (r) I t c i (estimated without radiation losses) have op- 
posite signs, the near equality of their magnitudes at 
t = r coo i in Figure 5 indicates that feedback and cooling 
are balanced within r coo i. A quasi-steady state is estab- 
lished with little or no net gas flow across r coo i , although 
feedback outflow continues in r > r coo i unaffected by ra- 
diation losses. 

The important agreement between opposing mass flows 
at r coo i shown in Figure 5 is a further confirmation of 
the self-consistency our simple feedback estimates. Less 
than 1 percent of the total feedback energy is stored as 
potential energy within the cooling radius. Nevertheless, 
central black holes - in their over-zealous, over-reaching 
efforts to feed back to the cluster atmosphere the energy 
acquired from gas cooling in their immediate vicinity - 
drive huge flows of gas out beyond the virial radius, but 
in the process necessarily provide enough mass outflow 
within the relatively small cooling radius to shut down 
the large cooling inflow that would otherwise occur in 
this critical central region. Although we do not consider 
the detailed time evolution of the initial adiabatic cluster 
gas profile as it transforms into the gas density profiles 
observed today, we imagine that this occurs in a quasi- 
steady manner, as explained above, in which feedback 
energy is widely distributed as PE in the cluster gas. 
By this means the feedback outflow always nearly bal- 
ances the radiative inflow, avoiding any large central gas 
concentration (and eventual overdensity due to star for- 
mation) or other excursions very far from the gas density 
profiles currently observed. 

2.9. Feedback production of cosmic rays 

It is likely that feedback consists of jets and jet- 
produced cavities that are filled mostly with cosmic rays. 
If so, it is interesting to compare the total feedback en- 
ergy in our composite clusters, ~ 10 63 ergs, with that 
expected from other cluster cosmic ray sources. For sim- 
plicity, consider proton cosmic rays that have cluster life- 
times comparable to the cluster age and compare the to- 
tal (proton) cosmic ray energy from feedback shocks to 
the much stronger accretion shock that produced the un- 
derlying entropy gradients in cluster atmospheres. From 
Table 2 the total potential energy of gas within the virial 
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radius in clusters 1 and 2 is \PE\ ~ 10 64 ergs so the to- 
tal thermal energy is E th ~ 0.5 x 10 64 ergs by the virial 
theorem. The virial temperature and entropy in the clus- 
ter gas are acquired in the accretion shock as the cluster 
formed. Typically, about 10% of the shock energy is con- 
verted to cosmic rays, so the total energy in cosmic rays 
created as the cluster formed is E cr ~ 0.1E t h ~ 0.5 x 10 63 
ergs. Therefore, the cosmic ray energy from feedback and 
cluster accretion are comparable. 

The thermal energy profiles in our adiabatic, pre- 
feedback cluster models (upper panels in Figs. 1-4) 
change very little after receiving the enormous feedback 
energies we consider. In our scenario almost all of the 
thermal energy created when feedback energy is initially 
deposited subsequently becomes potential energy as the 
cluster gas expands; this is the energy evolution de- 
scribed in Mathews & Brighenti (2008). More likely, 
most of the feedback energy is supplied from the central 
black holes as jets of supra-thermal cosmic rays that drive 
the shocks that increase the cluster gas entropy. When 
feedback is mostly in the form of cosmic rays, the net 
thermal energy of the cluster gas is subsequently reduced 
by cluster expansion more than it is initially increased 
by shock waves and the net effect of cavity production is 
to cool, not heat, the cluster gas (Mathews & Brighenti 
2008). AGN feedback energy and cavity production in 
cosmological calculations are assumed to be in the form 
of ultra-hot thermal gas. This injection of thermal en- 
ergy results in a small net increase in the total thermal 
energy of the cluster gas even after it has fully expanded. 
Since we do not explicitly consider the non-thermal cos- 
mic ray component here, a very small increase (rather 
than a decrease) in the post-feedback temperature pro- 
files appears in Figures 1-4. In either case, any small 
change in the thermal energy of the cluster gas can be 
neglected in our estimates of the total feedback energy. 

2.10. Locations of feedback energy deposition and 
storage 

It is important to recognize the distinction between the 
location in cluster gas where entropy-increasing feedback 
energy is deposited and where it is ultimately stored as 
potential energy which is what we address here. Evi- 
dently most feedback heating occurs in shock waves mov- 
ing away from advancing jets and expanding X-ray cav- 
ities. X-ray observations of these shocks generally indi- 
cate modest Mach numbers, Ai < 2, in which the gas 
entropy is increased only by < 1.2. The lower density of 
recently shock-heated gas results in an adiabatic, buoy- 
ant outward flow in the cluster that stops when the en- 
tropy of the heated region matches that of the ambient 
atmosphere. However, the entropy increases with clus- 
ter radius as S cx r mostly due to dissipation acquired 
during cosmic accretion. Consequently, gas heated with 
a M ~ 2 shock can only rise by a factor of ~ 1.2 in 
cluster radius where its feedback energy is stored as po- 
tential energy. Conversely, if the reduced gas fraction 
observed at 0.5r\, ~ 1 Mpc were due to shock heating 
events at 50 kpc (where X-ray cavities can be observed) , 
the entropy difference AS ~ 36 would require shocks 
with Mach numbers M ~ 17 in which the post-shock 
temperature would be increased by about 95. Shocks 
this strong have not been observed in cluster gas. Fur- 
thermore, it is likely that there is not enough cluster 



gas within 50 kpc which, when continuously heated this 
much with a filling factor of unity, can buoyantly supply 
the much larger mass of high entropy, feedback heated 
gas observed at 1 Mpc. In addition, mixing instabilities 
arc likely to dilute the entropy of heated gas during their 
long buoyant journey to r ~ 1 Mpc. Clearly, long range 
buoyant outflow is a losing proposition. The conclusion 
we draw from this is that feedback energy stored at some 
radius always exceeds the radius where the energy was 
initially deposited, but the difference between these two 
radii is unlikely to be very large. 

3. BLACK HOLE ACCRETION OR SPIN ENERGY? 

The time-averaged mechanical feedback power Lfb 
required to lift the adiabatic cluster atmospheres be- 
yond the virial radius r v can be estimated from Lfb ~ 
APE/td where t c i = 7 Gyrs is a typical age for 
large clusters. For composite cluster (1,2) the feed- 
back luminosity with and without initial cores is Lfi, ~ 
(0.38, 1.32) x 10 46 and L fb ~ (0.44, 1.50) x 10 46 erg s^ 1 
respectively. These luminosities are comparable to those 
of powerful quasars continuously active during the clus- 
ter lifetime. 

What are the implications of these enormous feedback 
powers for the mass of cluster-centered black holes? The 
conversion of accretion energy into black hole mass is 
more straightforward for luminous quasars in which the 
(thin) accretion disk can be directly observed. In this 
case, the total bolometric feedback energy radiated by 
an active black hole is related to the mass of the black 
hole by E ra d = e r MbhC 2 with e r ~ 0.1 (e.g. Davis & 
Laor 2010). While low redshift cluster-centered black 
holes generally have bolometric luminosities very much 
less than the estimated Lfb, it is unclear if this also ap- 
plies at higher redshift when most of the feedback energy 
may have been created. 

Suppose we adopt a mechanical feedback efficiency 
e m /i, such that the change in cluster gas potential en- 
ergy is related to mass accretion by 

APE = e mf bAM bh c 2 . (9) 

If we identify the total mass accreted with the mass of 
the black hole, AMbh ~ Mbh, and assume e m fb ~ 0.1, 
the resulting black hole masses for cluster 1, Mbh ~ 4.6 — 
5.4 x 10 9 Mq, are similar to those observed in cluster- 
centered galaxies, e.g. M bh w 6.4 x 10 9 M Q in M87 
at the center of the Virgo cluster (Gebhardt & Thomas 
2009). However, the corresponding black hole masses 
for the more massive cluster 2, Mbh ~ 1-5 — 1.6 x 10 10 
Mq, exceed those observed. This large black hole mass 
cannot be reduced by invoking the energy contributed 
to the cluster gas from accreting black holes in bulges 
of non-central satellite galaxies in the cluster. Almost 
all observed X-ray cavities and shock waves are associ- 
ated with cluster-centered black holes, not those in or- 
biting satellite galaxies. Moreover, the cluster center is 
where the virialized cluster gas is densest and has the 
shortest radiative cooling time. However, for cluster 
masses M v > 10 14 Mq, such as we consider here, the 
mass of the central galaxy (and therefore also its black 
hole) increases very little with increasing cluster virial 
mass M v (Lin & Mohr 2004). The mass of the central 
black hole in cluster A478 (in composite cluster 2) esti- 
mated from the bulge luminosity of the cluster-centered 
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galaxy, 5.8 x 10 9 M© (McNamara, Rohanizadegan, & 
Nulsen 2011), is considerably less than the mass esti- 
mated from M bh w APE(r v )/(e mfb c 2 ) « 1.6 x 10 10 M Q 
with e m fb — 0.1 for the cored cluster 2. Finally, the me- 
chanical accretion efficiency e m fb cannot be much larger 
than 0.1 and is likely to be much smaller. 

This difficulty can be alleviated if a substantial frac- 
tion of the required mechanical feedback energy is sup- 
plied by the rotational energy of central black holes 
which, when magnetically coupled to accretion disks, 
may form powerful feedback jets (Blandford & Znajaek 
1977). The maximum available energy from a rotating 
black hole, 0.29M bh c 2 = 5.3 x lO 62 [M b?l /(lO 9 M )] erg, 
is large enough to account for the gas depiction in both 
clusters 1 and 2. Recent models of the cosmological co- 
cvolution of galaxies and their black holes are consistent 
with the assumption that "radio mode" jet feedback de- 
rives from black hole spin and low-luminosity, advection- 
dominatcd central accretion (Sikora, Stawarz, & Lasota, 
2007; Fanidakis et al. 2011). There is considerable ev- 
idence for mechanical outflow from active galaxies and 
quasars (e.g. Crenshaw et al. 2003), but most relevant 
to our feedback estimates here are powerful FRII radio 
sources in luminous quasars that can transport enormous 
energies to great distances in the cluster gas (e.g. Mullin 
et al. 2008). Recently Singal et al. (2011) argue that the 
evolution of quasar luminosities at optical and radio fre- 
quencies are strongly correlated since redshift z w 3 with 
a significantly higher radio than optical evolution. In all 
likelihood, most of the radio-mode feedback energy from 
cluster-centered black holes probably occurred at earlier 
times. While spin energy is an attractive hypothesis, 
it must be reconciled with the approximately isotropic 
ejection of feedback jet energy currently observed in lo- 
cal galaxy clusters such as Virgo and Perseus and the 
rather large number of FRII radio sources with multiple 
hotspots, both of which, on the spin hypothesis, would 
indicate abrupt changes in the massive black hole spin 
axis on timescales that are implausibly small. 

Throughout this discussion we assume that the huge 
energy required to drive down cluster gas fractions orig- 
inates in (cluster-centered) black holes. While other as 
yet unidentified sources of energy cannot be excluded, a 
cluster-centered energy source seems likely. The largest 
bipolar X-ray cavities associated with feedback from clus- 
ter centers require energies ~ 10 62 ergs (McNamara et al. 
2005; Guo & Mathews 2010b) that are a significant frac- 
tion of the total feedback energies we estimate. Strong ad 
hoc preheating of baryons to a fixed adiabat before they 
flow across the cluster virial radius produces a character- 
istic entropy "floor" that may be inconsistent with the 
mass of gas observed in smaller group and galaxy po- 
tentials. To produce the quasi-powerlaw entropy profiles 
observed in clusters (dotted lines in Figs. 1-4) the pre- 
heated gas would need to be differentially heated in a 
fine-tuned manner, varying with redshift. 

4. RELATION BETWEEN GAS FRACTION AND ENTROPY 

In their recent paper Pratt et al. (2010) plot en- 
tropy and gas fraction profiles for over 30 nearby clus- 
ters observed with XMM-Newton, confirming the strong 
anti-correlation between these two parameters. When 
the mean entropy (S(r)) increases, the gas fraction (f g ) 
decreases and the amplitude of this behavior increases 



with decreasing cluster gas temperature and M v as ex- 
pected. But, in an interesting departure from previous 
work, in their Figure 9 Pratt et al. plot a combina- 
tion of the normalized entropy and the integrated, mass- 
weighted gas fraction [S(r)/ S 50 a](.f g (< r)/f b ) 2 / 3 with 
f/fsoo, showing that the previously discordant profiles 
for S (r) /^oo among their 30 clusters all collapse almost 
magically into a single tight bundle of correlations right 
along the relation expected for adiabatic baryonic col- 
lapse, S ad {r)/S 500 = 1.42(r/r 500 ) 1 - 1 for r < r 500 . 

It is easy to demonstrate a similar result with our 
composite clusters 1 and 2. The upper panels in Fig- 
ures 1 and 2 for "no core" clusters demonstrate that the 
gas temperature varies with cluster radius very much 
less than the observed or adiabatic entropy S(r) = 
T(r)/p(r) 2 / 3 in the central panels. The steady increase 
in entropy with cluster radius is due almost entirely to 
its reciprocal relation to decreasing gas density. More- 
over, the gas densities for the adiabatic and observed 
clusters vary with the total cluster density as p a d = fbpt 
and p bs = fgPt- Consequently, the adiabatic gas den- 
sity at any radius is related to the observed density 
by Pad = Pobs(fb/fg) where f b = 0.17. Using this 
relation, the local entropy in our adiabatic, feedback- 
free clusters can be approximately written in terms of 
the entropy in the observed cluster at the same ra- 
dius, S ad (r) = [T(r)/p(rf/\ bs (f g /f b f/\ where T(r) » 
T ad (r) » T obs (r). 

Dash-dotted lines in the central panels of Figures 1 
and 2 show profiles of S ad — S obs (f g /f b ) 2 ^ 3 which are 
almost identical to the solid line profiles for the adiabatic 
atmosphere, S a d(r). This is the result found by Pratt 
et al. Although we use the more appropriate local gas 

mass fraction f g (r) to show that S a d ~ S a d, the gas 
mass- weighted values of the integrated gas fractions f g (< 
r) used by Pratt et al. are not greatly different from 
local values f g (r) because the gas mass increases rather 
rapidly with cluster radius. Like Pratt et al. we find that 
the observed entropy profiles for both clusters 1 and 2 
collapse back to the adiabatic profiles when the observed 
density is increased by f b / f g . 

What is the physical significance of this result? It is not 
surprising that the gas densities between observed and 
adiabatic atmospheres are related by p a d{r) / p obs {r) = 
fb/ fg( r ) since this is required by the definition of these 
two atmospheres. However, the close agreement between 
S a d(r) and S a d(r) is possible because the gas tempera- 
ture profiles in the observed and adiabatic atmospheres 
(bound by the same virial mass) are virtually identical in 
clusters 1 and 2. To an excellent approximation, the ob- 
served and adiabatic entropy differ locally only because 
of the differing gas densities. It is also interesting that 
the entropy profiles of all clusters observed by Pratt et 
al. - and possibly all known clusters - can be modified in 
this way to recover the universal adiabatic profile. Clus- 
ters for which S a d(v) ~ S a d(r) are relaxed in the sense 
that non-gravitational "heating" events that increase the 
cluster gas entropy and drive it out in the cluster poten- 
tial, lowering the gas fraction, are either too weak or too 
old to retain evidence of increased gas temperature due 
to the most recent feedback event. Conversely, clusters 
for which S a d(r) slightly exceeds S a d(r) show evidence of 
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recent feedback-related increases in gas temperature and 
S a d(r) — S a d(r) is a measure of the location, energy and 
age of these feedback events. 

5. CONCLUSIONS 

By considering the difference between gas density pro- 
files in clusters created by non-radiating gravitational 
collapse and similar observed clusters, we estimate the 
total feedback energy received by the gas during the clus- 
ter lifetime. This estimate is insensitive to the precise 
time of the feedback events, but on average feedback en- 
ergy must be widely distributed in cluster radius. In- 
dividual feedback events produce jets and cavities and 
associated shock waves that heat the cluster gas. The 
heated, high entropy gas expands, ultimately leading to 
a expansion of the entire cluster atmosphere in which the 
feedback energy is stored as potential energy. The final 
cluster gas entropy profile is increased by feedback en- 
ergy and the gas fraction f g (r) profile is reduced below 
the cosmic baryonic value /& = 0.17. 

For clusters having masses M v > 10 14 M & the total es- 
timated feedback energy required to account for the ob- 
served depletion of cluster gas and the enhanced entropy 
profile is > 10 63 ergs, considerably in excess of cluster gas 
energy gained from supernovae or lost by radiation. Al- 
though enormous, this energy exceeds the largest known 
energy released in single, most powerful feedback events, 
< 10 ergs. The most likely source of this energy is 
feedback from central black holes in cluster-centered el- 
liptical galaxies. When averaged over a typical cluster 
lifetime, about 7 Gyrs, the mean mechanical luminos- 
ity, ~ 10 46 ergs s -1 , is comparable to that of powerful 
quasars. Such a large sustained luminosity may require 
energy creation not just from black hole accretion but 
also its spin energy. 

Immediately following a feedback energy event, we ex- 
pect a significant local increase in gas temperature and 
thermal energy. But after a cluster sound-crossing time 
the cluster expands and the temperature returns to a 
rather flat profile near the virial temperature required to 
support the cluster gas. The sound-crossing time to rsoo 
in our two composite clusters (1,2) is £500 = (4.6,3.5) x 
10 s yrs. In view of the insensitivity of the gas tempera- 
ture profile to the feedback energy after time ~ 10 8 yrs, 
changes in the cluster gas entropy S(r) — T(r)/p(r) 2 / 3 
arise almost entirely from changes in the gas density pro- 
file. After ~ 10 8 yrs, transient feedback events that 
increase the local gas temperature evolve by expansion 
into density reductions as thermal energy converts to po- 
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tential energy, retaining approximately the same global 
T(r). This explains the strong anti-correlation between 
excess entropy and reduced gas fraction in galaxy groups 
and clusters (Sun et al. 2009). Inspired by Pratt et al. 
(2010), we also demonstrate that the entropy profile ob- 
served in any relaxed cluster S(r), when multiplied by 
a factor containing the gas fraction f g (r), recovers the 
universal adiabatic gas entropy profile expected in the 
absence of feedback, S a d(r) = S(r)[f b / / g (r)] 2 / 3 cx r 1 ' 2 , 
with small deviations related to recent feedback events. 

Most of the feedback energy and entropy are delivered 
to very distant regions in cluster hot gas atmospheres, far 
beyond the cooling radius, where they have little or no 
effect on reducing the rate that gas cools near the central 
black hole, the presumed source of feedback energy. In 
the absence of feedback energy, the observed cluster gas 
is expected to cool at a rate M c f(^cooi) within the cool- 
ing radius r coo i, the cluster radius at which the cooling 
time equals the typical cluster age t c i. For our clusters 
fcooi ~ 100 kpc. Cooling flows near the central black hole 
can be greatly reduced or stopped if the mass of cluster 
gas that flows out across r coo i during the cluster lifetime 
t c i, driven by feedback, is approximately equal to the in- 
flowing mass due to radiation losses M c f{r coo i)t c i- For 
the clusters we consider it is gratifying that the cluster 
gas mass flow in both directions at r coo i - out due to 
feedback and in due to radiation losses - are very nearly 
equal, indicating that approximately the right amount of 
feedback energy is delivered to gas at r < r coo \ during t c i 
to drastically reduce cooling near the black hole where 
little or no cold or cooling gas is observed. Nevertheless, 
only a small fraction of the total feedback energy, less 
than 1 percent, is delivered and stored within r coo /. 

Finally, it is significant that the average mechanical 
feedback power Lfb ~ 10 46 erg s~ x implied by observed 
cluster gas fractions f g < fb — 0.17, are very substan- 
tially higher than estimates of Lft, from observations of 
X-ray cavities, ~ 10 59 ergs (Rafferty et al. 2006; McNa- 
mara, Rohanizadegan, & Nulsen, 2011). This discrep- 
ancy may be attributed to the difficulty of detecting X- 
ray cavities at distances exceeding about 50-70 kpc from 
cluster centers, particularly at higher redshifts, where 
most of the feedback energy is delivered and stored. 
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Fig. 1. — Two equilibrium atmospheres for composite cluster 1: a purely adiabatic virializcd gaseous atmosphere without central core 
(solid lines) and an observed atmosphere in a cluster of the same mass (dotted lines). Upper panel: Gas temperature profiles. Center 
panel: Gas entropy profiles. The dash-dotted line shows the adiabatic entropy 5 a£ j estimated from the observed entropy at the same radius 
(see Section 4). Lower panel: Gas density profiles. Points show observed gas densities for A133 and A383 from Vikhlinin et al (2006) 
which fit with the dotted line. The total cluster density profile pt(r) is shown with a dashed line. The cooling radius for cluster 1, 98 kpc, 
corresponds to logr coo [/r v i r = —1.30. 
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Fig. 2. — Identical to Fig. 1 but for composite cluster 2 with observed gas densities for A478 and A1413. The cooling radius for cluster 
2, 120 kpc, corresponds to \ogr coo i/r v i r = —1.35. 
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Fig. 3. — Two equilibrium atmospheres for composite cluster 1: a purely adiabatic virializcd gaseous atmosphere with central core (solid 
lines) and an observed atmosphere of the same mass (dotted lines). Panels arc identical to Fig. 1. 
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Fig. 4. — Identical to Fig. 3 but for composite cluster 2 with observed gas densities for A478 and A1413. 
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Fig. 5. — Sloping lines show profiles of the approximate time-averaged rate AM(r)/t c t that the local integrated gas mass is lowered due 
to feedback outflows at each (no core) cluster radius and i c ; = 7 Gyrs is the cluster lifetime. Horizonal lines show the rate M(r coo i) that 
gas cools by radiation as it flows inward across r coo ; in the absence of feedback. The cooling radius r coo ; is marked with vertical lines. 
(Dashed, solid) lines refer to composite cluster (1,2). For both clusters the mean rate of mass outflow at the cooling radius AM(r coo ;)/t c ; is 
very close to the mass inflow rate at this radius M c f( r cool) due to radiative losses. This near equality is a necessary condition for feedback 
energy to arrest cooling flows. 
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TABLE 1 

TWO COMPOSITE GALAXY CLUSTERS 3, 



cluster 


M 50 o b 


r$oo h 


T c 
± mg 


M v d 






(10 14 M ) 


(kpc) 


(keV) 


(10 14 M ) 


(kpc) 


A133 


3.17 


1007 


3.67 






A383 


3.06 


944 


4.37 






cluster 1 


3.11 


974 


4.02 


4.34 


1950 


A478 


7.68 


1337 


7.36 






A1413 


7.57 


1299 


6.81 






cluster 2 


7.63 


1319 


7.08 


11.3 


2682 



a Cluster data from Vikhlinin ct al. (2006). 



Cluster mass and radius at overdensity A — 500. 
c Density weighted mean temperature of cluster gas. 
d Cluster mass and radius at virial radius. 
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TABLE 2 
COMPOSITE GALAXY CLUSTER 
ATMOSPHERES 



r/r v M ? (r) a \PE\ A\PE\ 

(10 13 M ) (10 62 ergs) (10 62 ergs) 



CLUSTER 1: NO CORE 



ad b 


1.000 


7.382 


42.06 




ob b 


1.862 


7.383 


32.37 


9.69 


ad 


0.501 


4.459 


30.38 




ob 


0.933 


4.453 


24.17 


6.21 


ad 


0.269 


2.486 


19.57 




ob 


0.501 


2.476 


16.15 


3.42 



CLUSTER 2: NO CORE 



ad 


1.000 


19.22 


204.6 




ob 


1.549 


19.20 


171.6 


33.0 


ad 


0.501 


11.23 


142.4 




ob 


0.776 


11.19 


122.4 


20.0 


ad 


0.324 


7.362 


103.2 




ob 


0.501 


7.317 


90.26 


12.94 



CLUSTER 1: CORE 



ad 


1.000 


7.221 


40.46 




ob 


1.820 


7.275 


32.12 


8.34 


ad 


0.501 


4.299 


23.79 




ob 


0.891 


4.284 


23.58 


5.21 


ad 


0.282 


2.448 


18.74 




ob 


0.501 


2.476 


16.15 


2.59 



CLUSTER 2: CORE 



ad 


1.000 


19.00 


200.7 




ob 


1.549 


19.19 


171.6 


29.1 


ad 


0.501 


11.01 


138.5 




ob 


0.758 


10.96 


120.7 


17.8 


ad 


0.331 


7.325 


101.3 




ob 


0.501 


7.317 


90.26 


11.0 



a For each cluster/core combination, the integrated gas 
mass is chosen to be nearly the same for each pair of "ob" 
and "ad" solutions. 

"ad" and "ob" designate cluster atmospheres based on 
the adiabatic and observed gas density profiles. 



